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First-principles  study  on  Ni3Al  (111)  antiphase  boundary  with  Ti  and  Hf  impurities 
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The  effect  of  Ti  and  Hf  impurities  on  the  (1 1 1)  antiphase  boundary  (APB)  energy  of  Ni3Al  is  investigated  via 
ab  initio  calculations.  Cluster  expansion  is  performed  to  predict  supercell  total  energies  sampled  in  a  Monte  Carlo 
approach  that  accounts  for  nondilute  point  defects  at  finite  temperature,  obtaining  APB  energies  as  a  function 
of  impurity  concentration  and  temperature.  Of  the  two  ternary  elements,  Hf  is  more  effective  in  increasing  the 
APB  energy.  While  the  (111)  APB  energy  of  a  pure  Ll2  material  requires  at  least  second-nearest-neighbor 
interactions,  we  observe  a  strong  correlation  between  impurity-induced  APB  energy  enhancement  and  formation 
of  first-nearest-neighbor  Ni-Ni  bonds  across  the  APB  due  to  symmetry  breaking.  Using  a  linear-chain  model  and 
effective  bond  energies  derived  from  effective  cluster  interactions,  we  propose  a  mechanism  that  explains  why 
Hf  is  more  effective  than  Ti. 


DOI:  10.1 103/PhysRevB. 95. 214121 

I.  INTRODUCTION 

Ni3  A1  is  a  technologically  important  material  that  is  used  for 
y'  precipitate  strengthening  in  Ni-based  superalloys.  Further 
enhancement  in  strength  can  be  achieved  by  solute  additions, 
which  has  been  extensively  studied  in  experiments  [1].  Com¬ 
putational  techniques  for  solute-dislocation  interactions  have 
been  applied  successfully  to  predict  solute  strengthening  [2] 
and  solute  softening  [3].  In  precipitates  such  as  y'  significant 
strengthening  occurs  when  an  ordinary  (|[110])  dislocation 
cuts  the  precipitate  producing  an  antiphase  boundary  (APB). 
The  resistance  to  shear  is  proportional  to  the  APB  energy, 
which  is  sensitive  to  composition  and  temperature.  The  Ni-Al 
y'-y  phase  boundary  is  fairly  insensitive  to  temperature  and 
binary  composition,  so  the  ratio  of  Ni  to  A1  in  Ni3Al  is  fairly 
constant  over  a  range  of  alloy  compositions.  However,  ternary 
additions  can  strongly  affect  the  APB  energy  and  enhance 
strengthening  by  increasing  the  APB  energy.  Recently,  the 
effect  of  off-stoichiometry  [4]  and  that  of  transition-metal 
additions  [5,6]  on  the  APB  energy  of  Ni3Al  have  been 
investigated  from  computational  approaches. 

The  effect  of  substitutional  impurities  on  the  APB  energy  is 
a  subject  of  wide  interest.  Current  first-principles  approaches 
are  often  either  mean-field-based  (such  as  the  coherent 
potential  approximation  used  in  Refs.  [6-8])  or  atomistic  but 
not  statistical  (e.g.,  Refs.  [5,9]).  As  far  as  the  authors  are  aware, 
the  first  application  of  an  atomistic  and  statistical  method  to 
APBs  was  used  by  Asta  and  Quong  in  a  study  on  the  effect 
of  off- stoichiometry  in  TiAl  [10].  To  obtain  APB  energy  as 
a  function  of  temperature  and  impurity  concentration,  we 
have  developed  a  cluster  expansion  (CE)  method  based  on 
density  functional  theory  (DFT)  total-energy  calculations  and 
Monte  Carlo  (MC)  sampling  techniques  [11].  In  the  current 
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study,  we  apply  this  method  to  investigate  the  effect  of  Ti  and 
Hf  impurities  on  the  ^  ( 1 10)  {111}  APB  energy  in  Ni3Al  and 
propose  an  atomistic  mechanism  for  APB  energy  enhancement 
in  this  system. 

II.  METHODS 

DFT  [12,13]  calculations  within  the  Perdew-Burke- 
Ernzerhof  [14,15]  generalized  gradient  approximation  were 
performed  using  the  plane-wave  code  Vienna  Ab-initio  Sim¬ 
ulation  Package  (vasp)  [16-19]  with  projector  augmented 
wave  potentials  [20,21].  The  semicore  pseudopotential  was 
used  for  Hf.  To  achieve  accurate  phase  stability  in  the  Ni-Al 
system,  spin  polarization  was  applied  to  Ni  atoms  [22] .  Total 
energies  were  converged  to  within  10-6  and  10-4  eV  for  each 
self-consistent  loop  and  ionic  relaxation  step,  respectively, 
using  4000  Monkhorst-Pack  k  points  [23]  per  reciprocal  atom. 

For  the  cluster  expansion  (CE)  [24]  step,  71,  142,  and  107 
structures  were  used,  respectively,  for  the  binary  system  Ni-Al 
and  the  ternary  systems  Ni-Al-Ti  and  Ni-Al-Hf,  to  express 
the  total  energy  as  a  linear  combination  of  effective  cluster 
interactions  (ECIs): 

e= 53  •/“(n°r“)’  u) 

a  \  a  I 

where  the  summation  index  a  runs  over  clusters,  ma  and 
Ja  are  the  corresponding  multiplicity  and  ECI,  and  the 
bracketed  term  is  the  cluster  correlation  function.  A  6  x  6  x  5 
supercell  of  the  12-atom  (lll)-oriented  unit  cell  was  used 
for  the  MC  simulation,  performed  at  temperatures  between 
800  K  and  1600  K  and  impurity  concentrations  between  0.5% 
and  10%.  We  used  the  Alloy  Theoretic  Automated  Toolkit 
(atat)  software  package  [11,25,26]  for  CE  and  MC.  Further 
computational  details  can  be  found  in  Ref.  [11]. 

An  APB  has  the  desirable  property  of  not  breaking  the 
symmetry  of  the  system’s  underlying  lattice;  in  other  words, 
the  set  of  all  atomic  sites  (regardless  of  their  occupation)  forms 
the  parent  lattice  even  in  the  presence  of  the  APB.  As  a  result, 
a  standard,  translation-independent  cluster  expansion  can  be 
used  to  represent  the  energetics  of  the  system.  We  define  the 
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APB  energy  as 

E  -  2 E0  N(e  -  e0)  ... 

y  = - = - ,  (2) 

r  2  A  A 

where  eo  and  e  are  in  units  of  energy  per  atom,  N  is  the 
number  of  atoms,  and  A  is  the  cross-sectional  area  of  the 
supercell.  The  effectiveness  of  solute  additions  is  assessed  by 
the  enhancement  in  the  APB  energy: 

A y(x)  =  y(x)  -  y(0),  (3) 

where  x  is  the  impurity  concentration.  Contributions  from 
vibrational  entropy  have  been  found  to  be  small  in  this 
system  [4,27]  and  thus  are  not  considered. 

III.  RESULTS 

A.  Enhancement  of  APB  energy  by  impurities 

For  the  pure  M3AI  system,  its  equilibrium  lattice  constant 
from  DFT  and  APB  energy  from  CE  are  in  good  agreement 
with  the  literature,  as  reported  in  our  earlier  work  [11].  The 
ECIs  for  the  binary  Ni-Al  and  ternary  Ni-Al-Ti  and  Ni-Al-Hf 
systems  are  plotted  as  a  function  of  the  cluster  diameter  in 
Fig.  1.  The  CE  for  the  Ni-Al-Hf  system  has  a  cross-validation 
(CV)  score  of  33.0  meV/atom,  slightly  higher  than  that  for 
the  Ni-Al-Ti  system  (30.4  meV/atom).  The  APB  energy  of 
pure  N^Al  obtained  from  CE  of  the  binary  Ni-Al  system  is 
318  mJ/m2.  In  addition  to  the  agreement  with  the  unrelaxed 
first-principles  value  of  300  mJ/m2  reported  in  Ref.  [4]  versus 
344  mJ/m2  from  our  calculations,  we  validate  further  by 
relaxing  all  internal  degrees  of  freedom  at  constant  volume 
and  cell  shape.  We  obtain  an  APB  energy  of  217  mJ/m2, 


CD 

a 


pairs  triplets  quad. 


d/ dnn 


FIG.  1.  Effective  cluster  interaction  /  as  a  function  of  cluster 
diameter  d  normalized  by  nearest-neighbor  distance  dnn  for  (top) 
Ni-Al,  (middle)  Ni-Al-Ti,  and  (bottom)  Ni-Al-Hf  systems. 


I 

a 
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FIG.  2.  Effect  of  0.5%-10%  (left)  Ti  and  (right)  Hf  on  the  mean 
APB  energy  of  Ni3  A1  as  a  function  of  temperature. 


corresponding  to  an  error  of  16.4  meV/atom,  which  is 
consistent  with  the  CV  score. 

The  enhancement  of  APB  energy  by  impurity  X  [Eq.  (3)] 
is  obtained  from  CE  of  the  ternary  Ni-Al- A  system.  The 
reference  APB  energy,  y( 0),  is  248  and  259  mJ/m2  for 
Ni-Al-Ti  and  Ni-Al-Hf,  respectively,  which  has  a  much 
smaller  discrepancy  with  respect  to  the  DFT  APB  energy 
of  217  mJ/m2  than  the  binary  CE  reference  APB  energy  as 
reported  in  the  previous  paragraph.  To  formally  assess  the 
error  cancellation,  we  compare  the  difference  in  the  APB 
energy  of  pure  NpAl  obtained  from  CE  and  DFT  versus  that  of 
Ni3Al:Ti,  i.e.,  ()/“  -  yp^)  -  (y“  -  Values forthe 

first  bracketed  term  are  already  available:  y^re  =  248  mJ/m2 
and  yp^fj  =  2 1 7  mJ/m2 .  For  the  second  bracketed  term,  we  use 
a  2  x  2  x  2  (lll)-oriented  cell  with  composition  N^A^Ti, 
where  the  impurity  is  placed  near  the  interface.  The  APB 
energies  are  y^  =  287  mJ/m2  and  y^pT  =  253  mJ/m2. 
The  absolute  difference  between  the  two  bracketed  terms  is 
3  mJ/m2.  Although  this  value  may  not  reflect  the  accuracy 
of  CE  for  every  atomic  configuration,  it  is  indicative  of  the 
accuracy  of  our  predicted  effect  of  impurities  and  the  presence 
of  significant  error  cancellation.  Results  presented  henceforth 
are  based  on  ternary  CEs. 

Using  the  ternary  ECIs,  MC  sampling  of  the  mean  APB 
energy  change  is  plotted  in  Fig.  2  as  a  function  of  temperature 
for  each  impurity  concentration.  For  Ti,  maximum  enhance¬ 
ment  in  the  APB  energy  is  achieved  at  7.5  at.  %  within  the 
considered  temperature  range  (800  K  to  1600  K),  where  Ay 
is  above  50  mJ  m-2.  Adding  10  at.  %  Ti  is  less  beneficial  such 
that  Ay  is  even  smaller  compared  to  adding  5  at.  %  above 
1200  K.  On  the  other  hand,  Hf  is  predicted  to  enhance  the 
APB  energy  drastically.  A  concentration  of  2.5  at.  %  already 
outperforms  the  maximum  effect  of  Ti,  and  the  rate  of  increase 
in  APB  energy  gradually  reduces  as  Hf  concentration  increases 
to  10  at.  %,  where  the  enhancement  is  above  190  mJ  m-2.  This 
result  suggests  that  Hf  is  a  more  effective  candidate  than  Ti, 
which  agrees  with  experimental  findings  [28]. 
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In  the  following  subsections,  we  address  the  mechanism  by 
which  the  addition  of  Hf  is  so  effective  in  increasing  the  APB 
energy. 


B.  Correlation  to  impurity  coverage  near  APB 


While  the  simulations  presented  in  the  previous  section 
include  the  effects  of  all  possible  configurational  mechanisms 
on  the  APB  energy,  it  is  instructive  to  isolate  the  domi¬ 
nant  mechanism.  Atomistic  simulations  enable  a  formal  and 
detailed  investigation  of  this  question.  Our  general  approach 
is  to  identify  structural  features  that  best  correlate  with 
APB  energy  enhancement.  Specifically,  we  use  the  Pearson 
correlation,  defined  as 


cov(Ay,  X) 

PAy,X  =  - 

&Ay&X 


(4) 


where  o  is  the  standard  deviation,  to  further  investigate  the 
structural  effect  on  APB  energy  by  analyzing  the  correlation 
between  Ay  and  some  other  variable  of  interest,  X.  Here, 
values  of  p  range  from  —1  to  1,  where  p  =  0,  p  >  0,  and 
p  <  0  indicate  no  linear  correlation,  positive  linear  correlation, 
and  negative  linear  correlation,  respectively.  (Nonlinear  effects 
were  verified  to  be  negligible  by  fitting  to  a  fifth-order 
polynomial.)  The  Pearson  correlation  is  widely  used  in  the 
literature.  For  instance,  in  a  high-throughput  computational 
study  by  Pilania  et  al.  [29],  it  is  used  to  distinguish  strong 
versus  weak  correlation  between  various  physical  properties. 

In  this  subsection,  we  consider  the  coverage  of  impurity 
atoms  in  the  nearest  and  second-nearest  planes  to  the  APB, 
where  X  =  0\  and  02?  respectively.  Since  the  Pearson  coeffi¬ 
cient  is  invariant  under  scaling  and  constant  shift,  the  actual 
APB  energy  y  and  the  number  of  impurity  atoms  near  the  APB 
can  be  used  directly  in  Eq.  (4).  The  Ay  subscript  is  dropped 
from  p  hereafter. 

The  correlations  are  computed  from  all  the  sampled 
structures  in  MC  and  plotted  as  a  function  of  temperature  in 
Fig.  3  at  the  extreme  impurity  concentrations  (0.5%  and  10%). 


TABLE  I.  Correlation  between  APB  energy  and  the  coverage  of 
Ti  and  Hf  impurity  atoms  in  the  nearest  plane  (6\)  and  the  second- 
nearest  plane  (02). 


X  (%) 

T  (K) 

Ti 

Hf 

0i 

e2 

0i 

02 

0.5 

800 

0.63 

0.13 

0.75 

0.27 

1000 

0.62 

0.11 

0.73 

0.20 

1200 

0.60 

0.13 

0.72 

0.18 

1400 

0.53 

0.11 

0.68 

0.16 

1600 

0.40 

0.11 

0.60 

0.13 

1.0 

800 

0.63 

0.07 

0.77 

0.21 

1000 

0.62 

0.08 

0.77 

0.17 

1200 

0.63 

0.08 

0.74 

0.19 

1400 

0.57 

0.07 

0.69 

0.15 

1600 

0.48 

0.08 

0.66 

0.13 

2.5 

800 

0.60 

0.03 

0.78 

0.20 

1000 

0.61 

0.05 

0.77 

0.16 
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0.58 

0.04 

0.74 

0.17 

1400 

0.57 

0.05 

0.74 

0.14 

1600 

0.54 

0.06 

0.73 

0.14 

5.0 
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0.54 

-0.03 

0.67 

-0.01 
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0.53 

-0.07 

0.71 

0.06 
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0.53 

-0.06 

0.69 

0.04 

1400 

0.51 

0.02 

0.72 

0.05 

1600 

0.47 

0.01 

0.71 

0.01 
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0.47 

-0.10 

0.53 

-0.06 

1000 

0.42 

-0.09 

0.57 

-0.03 

1200 

0.43 

-0.05 

0.60 

-0.01 

1400 

0.41 

-0.05 

0.61 

-0.05 

1600 

0.39 

-0.07 

0.61 

0.03 

10 

800 

0.33 

-0.09 

0.31 

0.05 

1000 

0.31 

-0.07 

0.43 

-0.02 

1200 

0.31 

-0.07 

0.49 

-0.07 

1400 

0.27 

-0.06 

0.47 

0.00 

1600 

0.26 

-0.07 

0.54 

0.00 

Ti  Hf 


FIG.  3.  Correlation  between  APB  energy  and  the  coverage  of 
(left)  Ti  and  (right)  Hf  impurity  atoms  in  the  nearest  plane  (0i) 
and  second-nearest  plane  (02)  as  a  function  of  temperature.  For 
clarity,  only  the  correlations  at  0.5%  and  10%  impurity  concentrations 
are  shown.  Correlations  at  intermediate  concentrations  are  listed  in 
Table  I. 


Results  at  intermediate  concentrations  are  shown  in  Table  I. 
For  both  Ti  and  Hf,  coverage  of  the  nearest  layer  to  the  APB 
has  a  stronger  correlation  than  that  of  the  second-nearest  layer. 
The  correlation  value  pox  is  the  highest  at  0.5%  Ti  and  2.5% 
Hf  at  800  K,  being  0.63  and  0.78,  respectively.  In  general,  pex 
decreases  with  increasing  impurity  concentration  and  reaches 
a  minimum  at  10%.  This  result  suggests  that  impurity  coverage 
of  the  nearest  layer  is  only  somewhat  responsible  for  the 
enhancement  in  the  APB  energy,  especially  since  pex  can 
drop  to  values  around  0.3  at  10%  impurity  concentration,  and 
that  there  may  be  other  possible  mechanisms  for  APB  energy 
enhancement.  Correlation  between  the  enhancement  in  APB 
energy  and  impurity  coverage  of  the  second-nearest  layer  is 
weak,  since  \po2  \  is  within  0.3. 

We  remark  that  pex  is  not  predictive  since,  for  instance,  its 
peak  value  is  obtained  for  2.5%  Hf  at  800  K,  but  the  maximum 
Ay  occurs  at  10%  Hf  (Fig.  2). 

C.  Correlation  to  bonding  across  APB 

For  the  Ll2  structure,  the  number  of  each  type  of 
nearest-neighbor  bonds  is  invariant  upon  formation  of  the 
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FIG.  4.  Correlation  between  APB  energy  and  net  number  of 
first-nearest-neighbor  bonds  formed  across  the  APB  as  a  function 
of  temperature  for  (left)  Ti  and  (right)  Hf  at  0.5%  (open  symbols) 
and  10%  (solid  symbols)  concentrations.  Self-pairs  and  cross-pairs 
are  shown  in  the  top  and  bottom  panels,  respectively.  Correlations  at 
intermediate  Ti  and  Hf  concentrations  are  listed  in  Tables  II  and  III, 
respectively. 

|(ll0){lll}  APB.  Thus,  at  least  second-nearest  neighbors 
must  be  used  in  the  cluster  expansion  of  the  APB  structure 
of  pure  N^Al.  Here,  we  examine  the  net  number  of  first- 
and  second-nearest-neighbor  bonds  formed  across  the  APB. 
A  positive  (negative)  value  indicates  net  bond  forming 
(breaking)  between  two  species  X  and  Y  upon  formation  of 
the  APB.  Only  the  first-  and  second-nearest  planes  adjacent 
to  the  APB  are  considered.  In-plane  bonds  are  not  affected 
and  can  thus  be  ignored.  The  correlation  between  APB  energy 
and  first-  and  second-nearest-neighbor  bonds  across  the  APB 
is  shown  in  Figs.  4  and  5,  respectively,  as  a  function  of 
temperature  at  the  extreme  impurity  concentrations  (0.5%  and 
10%).  Results  at  intermediate  Ti  and  Hf  concentrations  are 
shown  in  Tables  II  and  III,  respectively. 

1.  First-nearest-neighbor  bonds 

While  first-nearest-neighbor  bonds  do  not  contribute  to  the 
(111)  APB  energy  for  pure  N^Al,  the  local  symmetry  in  the 
vicinity  of  the  APB  is  broken  when  impurities  are  introduced. 
Indeed,  from  Fig.  4,  the  magnitude  of  the  correlation  is  not  only 
different  from  0  but  also  very  high  in  several  cases.  We  examine 
all  possible  pairs  below  in  the  order  of  self-pairs  (Ni-Ni,  Al-Al, 
X-X)  followed  by  cross-pairs  (Ni-Al,  Al-X,  Ni-X). 

(a)  Ni-Ni.  For  both  NisAkTi  and  NisAkHf  systems,  strong 
correlation  is  observed  between  APB  energy  enhancement 
and  formation  of  Ni-Ni  first-nearest-neighbor  bonds.  In  the 
case  of  Ti,  the  correlation  is  weakly  dependent  on  impurity 
concentration  and  temperature,  ranging  from  0.72  (0.5% 


FIG.  5.  Correlation  between  APB  energy  and  net  number  of 
second-nearest-neighbor  bonds  formed  across  the  APB  as  a  function 
of  temperature  for  (left)  Ti  and  (right)  Hf  at  0.5%  (open  symbols) 
and  10%  (solid  symbols)  concentrations.  Self-pairs  and  cross-pairs 
are  shown  in  the  top  and  bottom  panels,  respectively.  Correlations  at 
intermediate  Ti  and  Hf  concentrations  are  listed  in  Tables  II  and  III, 
respectively. 

at  800  K)  to  0.91  (0.5%  at  1600  K).  In  the  case  of  Hf, 
the  correlation  becomes  stronger  with  increasing  impurity 
concentration,  reaching  a  maximum  of  0.87  (10%  at  800  K). 

(b)  Al-Al.  The  correlation  is  weak  in  both  systems,  where 
the  magnitude  is  within  about  0.4. 

(c) X-X.  In  the  case  of  Ti,  the  correlation  is  stronger  at  higher 
impurity  concentrations,  being  around  0.4  for  5%-10%.  In  the 
case  of  Hf,  it  is  the  highest  at  2.5%-5%,  being  around  0.5. 

(d)  Ni-Al.  For  Ti,  the  correlation  becomes  stronger  (more 
negative)  with  increasing  temperature  at  low  impurity  con¬ 
centration,  reaching  a  moderate  value  of  —0.665  (0.5%  at 
1600  K).  For  Hf,  the  correlation  generally  becomes  stronger 
with  decreasing  temperature,  reaching  a  moderate  value  of 
-0.667  (10%  at  800  K). 

(e)  Ni-X.  Moderately  strong  correlation  is  observed  for  the 
breaking  of  Ni-X  bonds  in  both  systems.  The  value  is  relatively 
insensitive  to  temperature  and  concentration,  being  around 
—0.6  for  Ti  and  —0.7  for  Hf. 

(f)  Al-X.  In  both  systems,  the  correlation  decreases  with 
increasing  concentration.  The  highest  value  is  0.62  for  Ti  (0.5% 
at  800  K)  and  0.71  for  Hf  (1.0%  at  800  K). 

For  all  the  first-nearest-neighbor  bonds,  the  correlation  is 
the  highest  for  Ni-Ni  bond  forming,  followed  by  Ni-X  and 
Al-X  bond  breaking. 

2.  Second-nearest-neighbor  bonds 

Compared  to  the  correlation  of  first-nearest-neighbor 
bonds,  that  of  second-nearest-neighbor  bonds,  in  general,  is 
opposite  in  sign  but  similar  in  magnitude. 
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TABLE  II.  Correlation  between  APB  energy  and  net  number  of  first-  and  second-nearest-neighbor  bonds  formed  across  the  APB  for 
Ni3Al:Ti. 


First-nearest-neighbor  bonds  Second-nearest-neighbor  bonds 


x(%) 

T(K> 

Ni-Ni 

Al-Al 

Ti-Ti 

Ni-Al 

Ni-Ti 

Al-Ti 

Ni-Ni 

Al-Al 

Ti-Ti 

Ni-Al 

Ni-Ti 

Al-Ti 

0.5 

800 

0.72 

-0.22 

0.21 

-0.08 

-0.63 

0.62 

-0.72 

0.21 

-0.25 

0.08 

0.64 

-0.63 

1000 

0.76 

-0.17 

0.17 

-0.14 

-0.62 

0.61 

-0.75 

0.17 

-0.21 

0.15 

0.62 

-0.62 

1200 

0.77 

-0.12 

0.20 

-0.20 

-0.60 

0.58 

-0.76 

0.11 

-0.22 

0.20 

0.60 

-0.59 

1400 

0.82 

0.05 

0.21 

-0.36 

-0.54 

0.52 

-0.81 

-0.07 

-0.21 

0.38 

0.54 

-0.53 

1600 

0.91 

0.41 

0.12 

-0.66 

-0.42 

0.41 

-0.89 

-0.44 

-0.13 

0.67 

0.43 

-0.42 

1.0 

800 

0.77 

-0.16 

0.28 

-0.15 

-0.63 

0.60 

-0.76 

0.17 

-0.29 

0.15 

0.64 

-0.63 

1000 

0.77 

-0.15 

0.27 

-0.15 

-0.62 

0.60 

-0.77 

0.15 

-0.30 

0.16 

0.63 

-0.62 

1200 

0.78 

-0.15 

0.25 

-0.18 

-0.64 

0.61 

-0.77 

0.14 

-0.30 

0.19 

0.64 

-0.62 

1400 

0.83 

-0.00 

0.23 

-0.33 

-0.58 

0.56 

-0.81 

-0.02 

-0.27 

0.35 

0.58 

-0.57 

1600 

0.87 

0.21 

0.21 

-0.51 

-0.51 

0.49 

-0.85 

-0.25 

-0.24 

0.53 

0.51 

-0.50 

2.5 

800 

0.81 

-0.07 

0.38 

-0.22 

-0.61 

0.54 

-0.80 

0.08 

-0.41 

0.23 

0.64 

-0.60 

1000 

0.80 

-0.08 

0.38 

-0.21 

-0.62 

0.54 

-0.80 

0.08 

-0.40 

0.23 

0.64 

-0.61 

1200 

0.82 

-0.03 

0.36 

-0.26 

-0.60 

0.52 

-0.81 

0.02 

-0.39 

0.29 

0.62 

-0.58 

1400 

0.83 

-0.01 

0.35 

-0.29 

-0.60 

0.53 

-0.82 

-0.01 

-0.38 

0.33 

0.62 

-0.59 

1600 

0.85 

0.08 

0.36 

-0.37 

-0.56 

0.50 

-0.84 

-0.11 

-0.38 

0.42 

0.59 

-0.55 

5.0 

800 

0.86 

0.12 

0.41 

-0.36 

-0.58 

0.41 

-0.86 

-0.13 

-0.46 

0.40 

0.63 

-0.55 

1000 

0.84 

0.07 

0.41 

-0.32 

-0.58 

0.43 

-0.84 

-0.09 

-0.44 

0.36 

0.62 

-0.56 

1200 

0.84 

0.07 

0.41 

-0.31 

-0.58 

0.43 

-0.84 

-0.09 

-0.44 

0.36 

0.61 

-0.54 

1400 

0.85 

0.09 

0.38 

-0.35 

-0.57 

0.43 

-0.85 

-0.12 

-0.41 

0.40 

0.60 

-0.55 

1600 

0.86 

0.13 

0.38 

-0.39 

-0.56 

0.44 

-0.85 

-0.18 

-0.39 

0.45 

0.59 

-0.55 

7.5 

800 

0.87 

0.20 

0.38 

-0.43 

-0.58 

0.38 

-0.88 

-0.26 

-0.44 

0.49 

0.63 

-0.55 

1000 

0.86 

0.22 

0.38 

-0.43 

-0.55 

0.35 

-0.87 

-0.29 

-0.41 

0.50 

0.59 

-0.49 

1200 

0.86 

0.16 

0.40 

-0.38 

-0.57 

0.39 

-0.86 

-0.24 

-0.43 

0.46 

0.60 

-0.50 

1400 

0.85 

0.15 

0.41 

-0.38 

-0.58 

0.39 

-0.86 

-0.23 

-0.41 

0.46 

0.61 

-0.53 

1600 

0.86 

0.18 

0.40 

-0.42 

-0.58 

0.40 

-0.86 

-0.25 

-0.41 

0.49 

0.62 

-0.55 

10 

800 

0.86 

0.26 

0.33 

-0.49 

-0.58 

0.37 

-0.87 

-0.40 

-0.40 

0.57 

0.62 

-0.52 

1000 

0.85 

0.24 

0.34 

-0.46 

-0.59 

0.38 

-0.87 

-0.38 

-0.39 

0.55 

0.61 

-0.51 

1200 

0.85 

0.25 

0.39 

-0.45 

-0.59 

0.33 

-0.86 

-0.36 

-0.41 

0.54 

0.61 

-0.49 

1400 

0.84 

0.19 

0.39 

-0.41 

-0.59 

0.36 

-0.84 

-0.31 

-0.41 

0.50 

0.61 

-0.50 

1600 

0.84 

0.20 

0.38 

-0.43 

-0.59 

0.39 

-0.84 

-0.31 

-0.42 

0.52 

0.62 

-0.53 

(a)  Ni-Ni.  Strong  correlation  is  observed  in  both  systems. 
For  Ti,  the  correlation  is  relatively  insensitive  to  concentration 
and  temperature,  ranging  from  —0.72  (0.5%  at  800  K)  to 
—0.89  (0.5%  at  1600  K).  It  is  more  sensitive  in  the  Hf  system; 
with  increasing  temperature,  it  becomes  weaker  (less  negative) 
at  high  concentration  but  stronger  (more  negative)  at  low 
concentration,  spanning  a  wider  ranger  from  —0.49  (1.0% 
at  800  K)  to  -0.89  (10%  at  800  K). 

(b) Al-Al.  For  Ti,  the  correlation  is  weak  since  its  magnitude 
is  within  about  0.4.  For  Hf,  it  is  weak  at  most  concentrations 
and  temperatures,  but  reaches  a  moderate  value  of  —0.69  at 
10%  at  800  K. 

(c)  X-X.  The  correlation  is  again  weak  for  Ti  with  a 
magnitude  of  about  0.4.  For  Hf,  it  is  moderate  at  2.5%-5%, 
approaching  —0.7. 

(d)  Ni-Al.  For  Ti,  the  correlation  is  moderate  (around  0.5) 
at  10%,  attaining  the  maximum  value  of  0.67  at  0.5%  at 
1600  K.  For  Hf,  the  correlation  is  the  strongest  (0.75)  at  high 
concentration  and  low  temperature. 

(e)  Ni-X.  For  both  systems,  the  correlation  is  relatively 
insensitive  to  temperature  and  concentration,  being  around  0.6 
for  Ti  and  0.7  for  Hf. 


(f)  Al-X.  Moderate  correlation  is  observed,  ranging  from 
about  —0.5  to  —0.6  for  Ti  and  about  —0.4  to  —0.7  for  Hf. 

Among  all  the  second-nearest-neighbor  bonds,  the 
strongest  correlation  is  found  for  Ni-Ni  bond  breaking.  All 
the  cross-pairs  (Ni-Al  forming,  Ni- X  forming,  Al- X  breaking) 
exhibit  moderate  correlation. 

D.  Mechanism  of  APB  energy  enhancement 
1.  Effective  bond  energy 

To  assess  the  atomistic  mechanism  of  APB  enhancement, 
here  we  define  an  effective  bond  energy  f  based  on  the  ECIs  of 
nearest-neighbor  pair  clusters.  We  show  how  it  can  be  applied 
in  the  next  subsection. 

From  Eq.  (1),  the  energy  of  a  cluster  a  is 


When  a  cluster  a  contains  only  a  pair  of  atoms,  m  and  n ,  each 
of  which  can  take  values  of  1  (Ni),  2  (Al),  or  3  (A),  we  use 
fmn  to  denote  the  effective  bond  energy. 
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TABLE  III.  Correlation  between  APB  energy  and  net  number  of  first-  and  second-nearest-neighbor  bonds  formed  across  the  APB  for 
Ni3Al:Hf. 


First-nearest-neighbor  bonds  Second-nearest-neighbor  bonds 


x(%) 

T  (K) 

Ni-Ni 

Al-Al 

Hf-Hf 

Ni-Al 

Ni-Hf 

Al-Hf 

Ni-Ni 

Al-Al 

Hf-Hf 

Ni-Al 

Ni-Hf 

Al-Hf 

0.5 

800 

0.51 

-0.41 

0.39 

0.16 

-0.75 

0.70 

-0.51 

0.40 

-0.51 

-0.16 

0.75 

-0.71 

1000 

0.54 

-0.39 

0.35 

0.13 

-0.74 

0.70 

-0.54 

0.38 

-0.45 

-0.13 

0.74 

-0.70 

1200 

0.56 

-0.34 

0.38 

0.09 

-0.73 

0.68 

-0.56 

0.33 

-0.48 

-0.07 

0.73 

-0.68 

1400 

0.63 

-0.21 

0.31 

-0.08 

-0.68 

0.65 

-0.62 

0.20 

-0.39 

0.10 

0.68 

-0.66 

1600 

0.73 

0.03 

0.27 

-0.33 

-0.62 

0.59 

-0.72 

-0.06 

-0.33 

0.35 

0.62 

-0.59 

1.0 

800 

0.49 

-0.43 

0.49 

0.20 

-0.77 

0.71 

-0.49 

0.42 

-0.60 

-0.19 

0.78 

-0.72 

1000 

0.56 

-0.40 

0.46 

0.15 

-0.77 

0.70 

-0.56 

0.39 

-0.57 

-0.13 

0.78 

-0.72 

1200 

0.56 

-0.36 

0.44 

0.11 

-0.75 

0.68 

-0.56 

0.35 

-0.55 

-0.09 

0.75 

-0.69 

1400 

0.59 

-0.26 

0.39 

0.00 

-0.70 

0.64 

-0.59 

0.24 

-0.49 

0.02 

0.71 

-0.66 

1600 

0.67 

-0.14 

0.36 

-0.15 

-0.68 

0.63 

-0.66 

0.10 

-0.46 

0.18 

0.68 

-0.64 

2.5 

800 

0.52 

-0.41 

0.51 

0.19 

-0.79 

0.68 

-0.50 

0.37 

-0.68 

-0.16 

0.80 

-0.70 

1000 

0.57 

-0.37 

0.53 

0.14 

-0.77 

0.65 

-0.56 

0.35 

-0.67 

-0.11 

0.79 

-0.69 

1200 

0.58 

-0.31 

0.51 

0.08 

-0.75 

0.62 

-0.58 

0.29 

-0.64 

-0.05 

0.76 

-0.66 

1400 

0.60 

-0.29 

0.53 

0.05 

-0.74 

0.63 

-0.60 

0.27 

-0.63 

-0.01 

0.77 

-0.67 

1600 

0.62 

-0.23 

0.53 

-0.02 

-0.74 

0.62 

-0.61 

0.20 

-0.63 

0.05 

0.76 

-0.65 

5.0 

800 

0.65 

-0.14 

0.43 

-0.09 

-0.70 

0.49 

-0.63 

0.02 

-0.64 

0.17 

0.74 

-0.52 

1000 

0.65 

-0.20 

0.50 

-0.03 

-0.74 

0.52 

-0.64 

0.12 

-0.68 

0.09 

0.77 

-0.54 

1200 

0.64 

-0.19 

0.49 

-0.03 

-0.72 

0.52 

-0.64 

0.12 

-0.66 

0.09 

0.75 

-0.55 

1400 

0.65 

-0.21 

0.55 

-0.01 

-0.75 

0.52 

-0.66 

0.14 

-0.68 

0.07 

0.77 

-0.58 

1600 

0.67 

-0.19 

0.54 

-0.05 

-0.74 

0.54 

-0.68 

0.11 

-0.67 

0.12 

0.77 

-0.59 

7.5 

800 

0.76 

0.04 

0.24 

-0.32 

-0.64 

0.44 

-0.76 

-0.28 

-0.55 

0.42 

0.67 

-0.41 

1000 

0.74 

0.01 

0.31 

-0.27 

-0.66 

0.45 

-0.74 

-0.20 

-0.57 

0.34 

0.70 

-0.44 

1200 

0.72 

-0.06 

0.39 

-0.18 

-0.69 

0.46 

-0.74 

-0.11 

-0.59 

0.28 

0.71 

-0.48 

1400 

0.72 

-0.05 

0.46 

-0.18 

-0.70 

0.44 

-0.74 

-0.10 

-0.63 

0.27 

0.72 

-0.46 

1600 

0.71 

-0.09 

0.49 

-0.13 

-0.71 

0.45 

-0.72 

-0.05 

-0.63 

0.23 

0.72 

-0.48 

10 

800 

0.87 

0.11 

-0.29 

-0.67 

-0.65 

0.69 

-0.89 

-0.69 

-0.32 

0.75 

0.64 

-0.57 

1000 

0.83 

0.12 

0.11 

-0.47 

-0.63 

0.54 

-0.84 

-0.48 

-0.44 

0.57 

0.63 

-0.47 

1200 

0.76 

0.01 

0.24 

-0.32 

-0.66 

0.50 

-0.79 

-0.33 

-0.50 

0.44 

0.65 

-0.43 

1400 

0.72 

-0.00 

0.32 

-0.27 

-0.65 

0.45 

-0.76 

-0.28 

-0.54 

0.39 

0.64 

-0.37 

1600 

0.70 

-0.04 

0.45 

-0.18 

-0.69 

0.41 

-0.73 

-0.18 

-0.59 

0.31 

0.68 

-0.39 

The  calculated  values  for  point  terms  and  first-nearest- 
neighbor  bonds  in  both  Ti  and  Hf  systems  are  given  in 
Table  IV.  We  emphasize  that  the  values  by  themselves  are 
not  meaningful  as  they  contain  point-term  ECIs  and  that  only 
differences,  obtained  by  computing  the  reaction  energy  at  a 
constant  number  of  each  atomic  species,  provide  desirable 
information.  In  the  next  subsection,  we  use  these  effective  bond 
energies  to  provide  a  quantitative  analysis  of  the  mechanism 
for  APB  energy  enhancement  and  the  differences  between  Ti 
and  Hf  impurities. 

2.  Mechanism 

In  Sec.  Ill  Cl,  we  have  shown  that  the  enhancement  in 
the  APB  energy  is  strongly  correlated  with  the  formation  of 
Ni-Ni  first-nearest- neighbor  bonds  across  the  APB.  This  result 


may  seem  unintuitive  at  first,  since  the  number  of  each  type  of 
first-nearest- neighbor  bonds  is  conserved  during  the  formation 
of  a  ^  ( 1 10)  {111}  APB  in  a  pure  Ll2  structure.  Moreover, 
consider  the  simple  linear-chain  reaction  model  with  only 
nearest-neighbor  interactions: 

Ni-Al-Ni  -*  Ni-Ni-Al,  (6) 

where  a  schematic  is  drawn  in  Fig.  6.  In  this  process,  the 
initial  state  has  an  energy  of  2/3i2  —  fc  and  the  final  state  has 
an  energy  of  fin  +  /fi2  —  /3\ .  The  reaction  energy  is  thus  f3\\  — 
Pn  ~  Pi  +  ft  -  Subsituting  the  values  from  Table  IV,  we  obtain 
—0.023  eV,  meaning  that  the  formation  of  Ni-Ni  first-nearest- 
neighbor  bonds  is  actually  energetically  favorable.  This  result 
is  in  agreement  with  the  finding  that  the  Ni-NisAl  interfacial 


TABLE  IV.  Calculated  ft  values  (in  eV)  for  point  terms  and  first-nearest-neighbor  bonds  in  the  Ni3Al: X  system. 


X 

ft  (Ni) 

ft  (Al) 

ft  (A) 

fti  (Ni-Ni) 

ft2  (Al-Al) 

ft3  (X-X) 

ft3  (Al-X) 

ft3  (Ni-X) 

ft2  (Ni-Al) 

Ti 

-0.481 

-  0.482 

-  0.062 

-  0.962 

-  0.926 

-  0.044 

-  0.604 

-0.563 

-  0.941 

Hf 

-  0.602 

-  0.425 

0.459 

-  1.143 

-  0.843 

1.063 

-0.012 

-  0.244 

-0.988 
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FIG.  6.  Schematic  of  the  reaction  in  Eqs.  (6)  and  (7)  projected 
onto  the  (1 1 1)  plane.  Solid  and  open  symbols  represent  atoms  in  the 
lower  and  upper  (111)  planes,  respectively.  Circles  (squares)  indicate 
Ni  (Al)  sites.  A  shaded  circle  (square)  in  the  upper  plane  represents 
an  impurity  X  on  the  Ni  (Al)  site,  whose  chain  is  drawn  in  thick  solid 
(dotted)  lines.  When  an  APB  is  created  along  the  [101]  direction, 
a  Ni-Ni  first-nearest-neighbor  bond  is  formed  in  both  XNi  and  XA1 
cases.  Equation  (6)  corresponds  to  the  change  in  the  dotted  chain 
without  impurity  substitution  on  the  Al  site. 


boundary,  across  which  Ni-Ni  bonds  are  formed,  can  have  a 
lower  energy  than  the  APB  energy  of  Ni3Al  [30]. 

However,  if  impurities  are  present  at  the  vicinity  of  the 
APB,  one  must  consider  another  possibility  (Fig.  6): 


TABLE  VI.  Reaction  energies  (in  meV)  for  site  substitution  by 
impurity  X. 


Xai  +  NiNi 

XNi  +  A1Ai 

X 

“ ^  XNi  +  NiAi 

Xai  +  AlNi 

Ti 

-250 

226 

Hf 

309 

-351 

obtained  as  follows.  Consider  the  substitutional  reaction 

X5+Aa^  Xa  +  A5,  (9) 

where  subscripts  A  and  B  represent  Ni  and  Al  sites,  respec¬ 
tively.  Note  that  in  the  Ll2  structure  the  A  site  has  4  B  and 
8  A  nearest  neighbors  and  that  the  B  site  has  12  A  nearest 
neighbors.  The  reaction  energy  is  expressed  conveniently  as 

A  Erxn  =  4^23  +  %Pi3  +  12/bi  —  (4^12  +  8/bi  +  12/3i3) 

=  4(/3n  +  ^23  —  P12  —  Pn)-  (10) 

(The  ECI  point  terms  cancel  out  exactly  in  this  procedure.) 
Equation  (9)  can  be  rewritten  more  abstractly  as 

(32)(11)  (31)(12).  (11) 


Ni-X-Ni  -*  Ni-Ni-X,  (7) 

with  reaction  energy  —  Pis  —  p>\  +  Ps,  which,  from  Ta¬ 
ble  IV,  is  equal  to  0.019  eV  for  Ti  and  0.162  eV  for  Hf. 
This  result  may  explain  why  the  addition  of  Hf  is  much  more 
effective  than  that  of  Ti  at  all  concentrations  (Fig.  2). 

On  the  other  hand,  the  formation  of  impurity  X-X  bonds 
across  the  APB  may  seem  to  be  a  more  intuitive  mechanism 
in  the  sense  that  the  overlap  of  impurity-induced  elastic 
fields  should  be  energetically  unfavorable.  This  is  verified  by 
considering  the  linear-chain  model: 


Compared  with  Eq.  (10),  the  reaction  energy  is  proportional 
to  the  difference  between  the  sum  of  inner  and  outer  pairs 
(right-hand  side:  11,  23;  left-hand  side:  12,  13).  Similarly,  the 
other  substitutional  reaction, 


XA  +  Bb  —>  Xb  +  Ba, 
(31)(22)  ->  (32X21), 
has  a  reaction  energy  of 

A£rxn  =  4(^22  +  P\3  —  Pl2  —  fe)- 


(12) 


(13) 


X-H-X  ->  X-X-H,  (8) 

where  H  represents  a  host  atom  (Ni,  Al).  The  reaction  energy 
is  Pss  —  P13  —  Ps  +  Pi  and  is  summarized  in  Table  V,  which 
shows  that  Hf  atoms  tend  to  avoid  each  other  more  than 
Ti  atoms  do.  While  the  formation  of  the  Hf-Hf  bond  is 
indeed  more  energetically  unfavorable  than  that  of  Ti-Ti,  its 
correlation  to  APB  energy  enhancement  is  only  moderate  when 
compared  to  that  of  Ni-Ni  (Fig.  4). 


The  correspondence  between  the  index  pairs  and  substitutional 
reaction  energy  is  a  special  property  of  the  Ll2  structure  due 
to  the  ratio  of  nearest-neighbor  species. 

Using  the  effective  bond  energies  in  Table  IV,  the  site 
preference  of  Ti  and  Hf  impurities  in  Ni3Al  are  obtained 
from  the  above  energetics  analysis.  Reaction  energies  for  site 
substitution  are  given  in  Table  VI.  By  comparing  the  sign 
of  reaction  energies,  the  Ni  (Al)  site  is  preferred  for  Ti  (Hf) 
substitution. 


E.  Site  preference 

As  another  application  of  effective  bond  energies,  the  site 
preference  of  an  impurity  X  in  the  Ll2  host  A3B  can  be 


TABLE  V.  Reaction  energies  (in  meV)  for  formation  of  X-X 
first-nearest-neighbor  bonds  in  Eq.  (8). 


X 

Host  atom 

Ti 

Hf 

Ni 

100 

246 

Al 

140 

191 

IV.  DISCUSSION 

This  study  involves  two  main  parts:  (1)  modeling  of 
impurities  and  (2)  investigating  the  mechanism  of  APB  energy 
enhancement.  In  this  section,  we  discuss  the  highlights  and 
limitations  of  our  approach  from  these  two  perspectives. 

Regarding  the  modeling  of  impurities  for  APB  energetics, 
typical  methods  use  total-energy  calculations  of  a  few  selected 
structures  (e.g.,  Ref.  [9])  or  the  coherent  potential  approxima¬ 
tion  (e.g.,  Ref.  [6]).  Due  to  the  vast  number  of  possible  struc¬ 
tures,  the  former  method  cannot  provide  accurate  statistics 
with  current  computing  capabilities.  The  latter  method,  being 
a  mean-field-type  technique,  avoids  the  statistical  problem, 
but  trades  off  with  atomistic  details.  In  contrast,  the  CE-MC 


214121-7 

Distribution  A.  Approved  for  public  release  (PA):  distribution  unlimited. 


SUN,  WOODWARD,  AND  VAN  DE  WALLE 

approach  used  in  this  study  is  both  statistical  and  atomistic  by 
definition,  but  the  drawback  is  that  long-range  elastic  effects 
cannot  be  captured.  Nonetheless,  the  APB  energy  predicted 
from  CE  is  an  important  quantity  that  can  be  validated  by 
comparing  to  total-energy  calculations  [11]. 

As  for  the  mechanism  of  APB  energy  enhancement,  the 
coverage  of  impurities  near  the  APB  is  naturally  assumed 
to  be  responsible  for  such  an  effect.  Indeed,  if  the  impurity 
were  located  deep  within  the  bulk,  the  APB  energy  would 
not  be  affected.  Given  that  first-nearest-neighbor  bonds  are 
not  involved  in  (111)  APBs  of  pure  Ll2,  impurity  coverage 
at  the  APB  is  necessary  for  APB  energy  enhancement,  but 
the  exact  mechanism  requires  deeper  investigation.  Using  a 
statistical  approach,  correlations  with  various  quantities  can  be 
examined  in  detail.  In  particular,  Ni-Ni  first-nearest-neighbor 
bond  formation  across  the  APB  has  a  higher  correlation  to  APB 
energy  enhancement  than  impurity  coverage  does,  revealing 
that  local  structure  (i.e.,  arrangement  of  atoms)  can  play  a  more 
important  role  than  local  chemistry  alone  (i.e.,  atomic  species). 
The  mechanism  proposed  in  Eq.  (7)  is  a  result  of  comparing 
all  first-nearest-neighbor  bonds.  While  certain  second-nearest- 
neighbor  bonds  also  have  high  correlation  to  APB  energy 
enhancement,  their  ECIs  are  much  smaller  in  magnitude 
(Fig.  1)  than  first-nearest-neighbor  ECIs.  Triplet  ECIs  in  both 
ternary  systems  are  comparable  in  magnitude  to  first-nearest- 
neighbor  pair  ECIs,  suggesting  that  local  structures  of  higher 
complexity  may  influence  the  APB  energy  as  well.  To  tackle 
the  general  problem  of  identifying  complex  structures  at  the 
vicinity  of  the  APB,  a  machine-learning  approach  is  suitable 
for  further  study.  Although  simplistic,  the  linear-chain  model 
of  reaction  energies  presented  in  Sec.  HID  2  is  still  useful  for 
a  first- step  analysis  of  other  materials. 

The  formal  procedure  of  using  effective  bond  energy  to 
predict  site  preference  has  been  outlined  in  Secs.  IIID1 
and  III  E.  We  do  not  attempt  to  resolve  the  conflicting  results 
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of  site  preference  in  the  literature  (e.g.,  reported  in  Ref.  [1]) 
as  some  of  these  differences  may  be  the  result  of  considering 
different  compositions  of  the  host  structure. 

Finally,  we  note  that  theoretical  studies  of  mechanical  be¬ 
haviors  such  as  yield-stress  anomaly  and  ductility  often  involve 
ratios  of  APB  energies  corresponding  to  various  slip  systems, 
which  are  used  in  stability  maps  derived  from  mesoscale 
dislocation  theory  [31,32].  The  method  presented  in  this  paper 
can  be  extended  to  investigate  the  temperature  dependence  or 
the  concentration  dependence  of  these  phenomena  from  first 
principles. 

V.  CONCLUSIONS 

We  have  presented  a  first-principles  Monte  Carlo  approach 
to  investigate  the  effects  of  Ti  and  Hf  impurities  on  APB 
energies  of  N^Al.  The  (111)  APB  energy  is  increased  by 
both  solute  additions,  where  Hf  addition  leads  to  stronger  en¬ 
hancement.  APB  energy  enhancement  is  positively  correlated 
to  impurity  coverage  of  the  APB,  and  a  strong  correlation  to 
formation  of  first-nearest-neighbor  Ni-Ni  bonds  is  observed. 
Based  on  a  linear-chain  model,  the  enhancement  mechanism 
is  proposed  as  an  impurity-activated  process,  Ni-X-Ni  -> 
Ni-Ni-X,  where  the  reaction  energy  is  19  meV  for  Ti  and 
162  meV  for  Hf. 
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